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Abstract. High momentum jets and hadrons are important probes for the quark gluon plasma 
(QGP) formed in nuclear collisions at high energies. We investigate how fluctuations in the 
background density of the QGP and fluctuations in the spatial distribution of the hard process 
create azimuthal asymmetries of the high momentum hadron spectrum, described by the Fourier 
coefficients Vn, n > . We estimate the coefficients up to va in a simple energy loss model tuned 
to single inclusive hadron suppression. 



With the study of quark gluon plasma in collisions of relativistic heavy ions moving into 
an increasingly accurate quantitative phase it has been found important to include fluctuations 
in the space-time structure of the fireball into calculations of bulk quantities [H [2]. These 
can emerge from fluctuations in the initial energy density e{x,y) in the plane transverse to 
the beam axis and leave signature effects on bulk quantities like the azimuthal asymmetry 
coefficients Vn- For example fluctuations can lead to a sizeable triangular flow vs which would 
be vanishing in an averaged fireball due to the overall geometry of the nuclear overlap [2]. At 
large transverse momentum Px fluctuations in the position of the hard process can also affect 
observables accessible in current heavy ion experiments. 

Here we have explored the role of such fluctuations on the suppression of high-Pj^ hadrons 
and their generalized azimuthal asymmetry coefficients f„. The Vn are defined as a Fourier 
decomposition of the azimuthal angle dependent spectrum 



where the angle $ is measured with respect to the reaction plane defined below and the 5„ 
are phases that encode a misalignment with the reaction plane. We note that for smooth, 
non- fluctuating fireballs we expect all odd coefficients vi, etc. at midrapidity to vanish for 
symmetry reasons. 

However, in any given single event the initial energy density will typically exhibit a non- 
vanishing triangular eccentricity €3 which could in turn lead to a non- vanishing V3 . The event-by- 
event fluctuations in initial energy density are driven by fluctuations of the positions of nucleons 
in the initial nuclei and in the amount of energy deposited around midrapidity for every nucleon- 
nucleon collision [3]. We also expect that the triangular eccentricity is not correlated with the 
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Figure 1. Engineered events witli n = 3, 4, 5, 6 

reaction plane, i.e. 83 should appear random. Similar arguments can be made for other n > 
and we expect all odd Vn to acquire non-vanishing values with realistic fluctuations. While all 
of this has first been discussed for the bulk of the fireball [3j these statements easily transfer to 
hard probes. Fluctuations in the energy density lead to fluctuations in energy loss. In addition, 
the position of a hard process which creates a hard probe is subject to fluctuations. 

Systematic measurements of w„ at large momentum could lead to further constraints on the 
type of energy loss prevalent in QGP, and on the size of the transport coefficient q. It can also 
give an independent handle on the size and granularity of initial state fluctuations. Here we 
report on a quantitative study of high momentum azimuthal coefficients using a simple energy 
loss model with realistic fluctuations. 

We calculate high momentum hadron spectra using our simulation package PPM [H [5]. It 
samples initial momentum distributions of quark and gluon jets from a perturbative calculation 
and propagates leading partons through a given background fireball. Different energy loss models 
can be employed. Here we will show results from a simple LPM-inspired (sLPM) deterministic 
energy loss model dE /dx ~ qx where q scales with the 3 /4th power of the local energy density [4j . 
As a cross check we will sometimes also use the non-deterministic Armesto-Salgado- Wiedemann 
(ASW) model [61 [7j. Fitted to the same experimental data on single hadron suppression these 
models cover a wide range of values for q. PPM will eventually fragment leading partons into 
hadrons and all results here will be shown for pions. 

We have used the Glauber Monte Carlo generator GLISSANDO [8] to produce an ensemble 
of Au+Au events at top RHIC energy using the three different centralities b = 3.2, 7.4 and 
11 fm. However, first we check the general relation between the spatial eccentricities [9] 
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and the azimuthal asymmetry coefficients Vm using "engineered" events with particular fixed 
eccentricities. These are created with the energy density modeled as simple Gaussians in the 
transverse plane with a cos n(p undulation of the mean square radius, as shown in the examples 
in Figfl} 

First we scanned the space (fni £m)) TijiTi — 1, 6 in search of correlations. As expected we find 
non-zero Vn for a given only if n — ifn where i > is an integer. Next we explored the scaling 
of Vn with the size of the eccentricity £„. We expected a monotonically increasing function Vn{en) 
which can be seen confirmed in Fig. [2} Deviations from monotony only occur for unrealistically 
large eccentricities. These basic results should hold if realistic fluctuations are considered. Fig. [3] 
shows the correlation between V2 and 62 for our ensemble of GLISSANDO events for all 3 impact 
parameters for Au-|-Au collisions at top RHIC energies. The basic linear correlation persists for 
n = 2 but is washed out. Correlations for n > 2, not shown here, are much more weakened to a 
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Figure 2. Azimuthal asymmetry Vn vs eccentricity e„ in engineered events for n = 3, 4, 5, 6. 
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Figure 3. Correlation between V2 and €2 for an ensemble of Au+Au collisions at top RHIC 
energy in GLISSANDO for three impact parameters 3.2, 7.4 and 11 fm. 



point that makes it hard to predict Vn for a known e„. All results shown are for sLPM energy 
loss but the corresponding result for ASW show no noticeable difference, making the conclusions 
rather robust against large variations in the microscopic origin of energy loss. 

We determined the phases 5n for high momentum pions in our ensemble of GLISSANDO 
events for three different impact parameters. In Fig.|4] we plot those phases relative to the 
reaction plane determined by the eccentricity €2- In other words our definition of a reaction 
plane is given by the fundamental initial ellipticity, which in general deviates from the plane 
defined by beam axis and impact vector. We observe that 63 and 5^ are randomly distributed, 
so there is no correlation between the reaction plane and the fiuctuations which create €3 or £5. 
This had also been found for the bulk azimuthal asymmetries before [9j . ^4 and 6q on the other 
hand show a correlation with the reaction plane which is consistent with V4 and vq receiving 
contributions from £2- 

The transverse momentum dependence of the coefficients Vn for pions for our ensemble of 
GLISSANDO Au+Au events is shown in Fig. [5j We observe two hierarchies of coefficients, one 
being V2 > V4 > vq and the second one being vi > V3 > v^. V2 is always the largest coefficient, 
even in the most central events. Interestingly vi is non-zero and the second largest coefficient, 
beating and V4 by more than a factor 2. Momentum conservation dictates a sum rule for vi 
integrated over Pt- The recoil of the medium in which energy is lost would lead to a negative vi 
at lower momentum. Such a back reaction is not included in this calculation, vi at intermediate 
and large Pt could be sensitive to the mechanism of medium recoil. Generally we point out 
that our results start to become unreliable below 4 to 6 GeV since hydrodynamic expansion 
was not included in the calculation. Fig. [5] shows the results for both sLPM and ASW energy 
loss. We have to conclude that azimuthal asymmetry coefficients are not particularly useful to 




Figure 4. Distribution of the phases 5n for n = 3, 4, 5, 6 in our ensemble of Au+Au events. 
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Figure 5. Pion Vn vs Pt for two impact parameters for Au+Au colhsions at top RHIC energy, 
calculated with either sLPM or ASW energy loss. 

discriminate between energy loss models. We also find that the Pt dependence of the coefficients 
becomes rather weak at large momenta. 

To summarize, in this study we have explored higher order azimuthal asymmetry coefficients 
at large momentum in heavy ion collisions. We find that in general f„ rises with but this 
correlation weakens for larger n. We also find that there are only a few cross correlations 
between Vn and for n ^ m. We have also classified the preferred angular orientation of the 
azimuthal asymmetries as given by the phases 5„ with respect to the reaction plane. We find a 
decorrelation with the reaction plane for all odd n. Finally we have made predictions for Vn as 
a function of Pt in two different energy loss models. We find mostly consistent results between 
those two models with V2 being the largest coefficient followed by vi. In general the Vn carry 
geometrical information and measurements at large momentum can be complementary to those 
for bulk observables, but they seem less useful to distinguish different energy loss models. 
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